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Summary 


A pulsed neutron source was introduced into the Swedish heavy-water reactor Rl when it 
was in a subcritical state. The time behaviour of the neutron density was studied by means of a 
multichannel time analyser. Using this technique the time constant of the prompt neutrons and 
the reactivity was determined at several heavy water levels down to about 20 cm below the critical 
level. It was found that the reactivity equivalent of the heavy water is (7.2+0.1) x 10-5 per mm, 
and that the mean lifetime of the thermal neutrons in the reactor is 0.68+0.02 msec. These 
values refer to the critical level and are based on a fraction of delayed and photo neutrons amount- 
ing to 0.78 x 10~?. A calibration of the safety rods gave the result that each of them corresponds 
to a reactivity of 1.7+0.2 per cent and both together to 3.4+0.4 per cent. 


Introduction 


Among the quantities which are of fundamental importance to know for the 
operation of a reactor are the reactivity equivalents of the safety and control rods. 
In a heavy-water reactor the dependence of the reactivity on the amount of heavy 
water is also of great importance. 

For the determination of reactivities of this kind two methods are generally used. 
The first method is to measure the period of the supercritical reactor and calculate 
the excess reactivity with Nordheim’s formula [1, 2, 3]. In the second method the 
influence of different absorbers on the multiplication of the subcritical reactor is 
studied by measuring the neutron density with a counter [3, 4]. The advantage is 
that large reactivities can be measured without the hazards involved in having the 
reactor supercritical with a large excess reactivity, but the method gives only relative 
values. In order to obtain the absolute reactivities at least one measurement according 
to the first method must be made. 

The mean lifetime of the thermal neutrons in a reactor is of interest for the kinetic 
behaviour and for estimating the source strength of the reactor. It can be determined 
by oscillating an absorber in the reactor and measuring the phase and amplitude of 


1 A preliminary report on this work was given at the Geneva Conference, paper P/789 (1955). 
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the variations caused in the neutron density as a function of the frequency [5, 6, 7, 
8]. Another method is to investigate the time behaviour of the neutron density after 
a step change in the reactivity [8, 9]. 

There is, however, a method (suggested by G. von Dardel) by which it is possible 
to determine both reactivities and the neutron lifetime from measurements on a 
subcritical reactor. In this method the reactor is pulsed with a neutron source and 
the time behaviour of the neutron density is studied. In the present paper such experi- 
ments with the Swedish heavy-water reactor R1 will be described. Similar experi- 
ments have independently been made in Russia on uranium-graphite lattices [10]. 


Theory 


Suppose, that a source inside a subcritical reactor produces rectangular neutron 
pulses, each containing S neutrons. It is assumed for simplification of the treatment 
that the repetition time is much shorter than the shortest periods of the delayed 
neutrons and that the time constant of the reactor is much shorter than the pulse 
length. The neutron density in the reactor will then vary with time as shown in Fig. 1. 
The prompt multiplication gives rise to pulses of the same shape as the source pulses. 
In addition there is a constant density due to the delayed neutrons, which cannot 
follow the rapid variations of the source. The area of a pulse above the constant 
value in Fig. 1 is denoted by A, and the area under the delayed neutron curve during 
one period by A,. The prompt multiplication and the total multiplication respectively 
correspond to [11] 


a:S 
A sok 1 
TSE) (1) 
a:S 
Ad te 2 
eit: 1 — kere 2) 


where f is the fraction of delayed neutrons (for a heavy-water reactor including photo 
neutrons) and a is a constant. Equations (1) and (2) give directly 


A, ¥. 1 — kere ev. 
a, phe oe (3) 


< 


where 9, the negative reactivity, is defined as 


= Le Kerr 
: Kese (4) 


If the fraction of delayed neutrons is known, it will thus be possible to determine 
the reactivity by a measurement with a pulsed neutron source. 

In the case actually considered the simplifying assumptions made above are not 
valid. The finite value of the reactor time constant, | p and the short periods of delayed 
neutrons modify the curve to a shape shown in Fig. 2. However, the area under the 
curve can still be divided into two parts, one arising from the prompt and one from 
the delayed neutrons. It is shown in Appendix I that with the usual approximations 
for a bare thermal reactor [11], assuming the age of the delayed neutrons to be the 
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Fig. 1. The neutron density in a pulsed sub- Fig. 2. A measured curve obtained with the 
critical reactor as a function of time under reactor subcritical by about one per cent. 
the simplifying conditions given in the text. The pulse length was 0.25 sec and the repe- 


tition time 0.5 sec. Measuring time was one 

hour. The dotted line indicates the neutron 

density arising from delayed neutrons. Ac- 
cording to eq. (3), A,/A, =0/f. 


same as the age of the prompt neutrons, the ratio of areas is still given by eq. (3). 
The two-group theory for a reactor with reflector also yields the same result. 

The rise and decay time of the neutron density is determined by the reactor time 
constant /,,. From eq. (Al) in Appendix I we see that 


f l _— U(1+e) 
POATEAT= Ete Bi () 


where / denotes the mean lifetime of the thermal neutrons in the reactor. By measuring 
the time constant and using the reactivity value obtained from eq. (3) the mean 
lifetime can thus be determined. 


Experimental arrangement 


The measurements were made with the Swedish heavy-water reactor Rl in a 
subcritical state [12]. A schematic view of the experiment is seen in Fig. 3. The 
neutron pulses were produced by means of a 150 kV acceleration tube [13]. It was 
placed horizontally in front of the reactor shield, and a 4 m long aluminium tube, 
5 cm in diameter, was extended into the reactor tank through a horizontal channel. 
A zirconium-—tritium target was attached to the end of the tube. Deuterium ions 
striking the target produced 14 MeV neutrons, by the T—D reaction. By modulation 
of the ion source [14] a pulsed neutron source was thus obtained about 30 em from 
the centre of the reactor. 

For convenient control of the ion beam it was passed through a hole in a mirror 
at 45° to the beam direction. The beam was strongly defocussed:during its passage 
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to the target, because the pumping velocity in the target tube was rather low. 
The tube walls received about as large an ion current as the target and this gave 
rise to unwanted D-D neutron sources. However, the broadness of the beam had 
the advantage that higher current could be allowed to reach the delicate target 
material without any risk of damaging it. 

In another horizontal channel which is about 20 cm lower than the target channel, 
a small BF, proportional counter mounted on a 4 m probe was introduced. The 
sensitive volume was 2.5 cm? and the counter was filled with 50 mm BF, and argon 
to a total pressure of 700 mm Hg. After amplification the pulses from the counter 
were fed into a 40-channel time analyser [15], by which the time behaviour of the 
neutron density in the reactor could be studied. 

Channel widths of the order of 0.01 sec were used for the measurements. They 
were obtained by feeding the pulses from the second frequency divider through an 
external scaler to the gate generators. The curve shown in Fig. 2, for example, was 
taken with a channel width of 4 x 12 x 256 = 12288 usec. 

The first 20 channels of the time analyser feed scale-of-16 units and the other 20 
feed scale-of-4 units. With the present channel widths the counting rate which could 
be used was limited by the mechanical registers connected to the scalers. To avoid 
counting losses in the registers the time analyser was coupled so that only every 
second channel was used, but with a scale-of-64 reduction. 


Reactivity measurements 


For the determination of the reactivity equivalent of the heavy water, measure- 
ments were made at 7 different water levels down to 18 cm below the critical level. 
All safety and control rods were out of the reactor. The pulse repetition time was 
0.5 or 1 sec and the pulse length half this period. The target current during the pulse 
was between 10 and 20 wA. The neutron production was always started at least 10 
min before the measurement commenced in order to reach equilibrium as regards 
the delayed neutrons. 
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Table 1. 

H & UB 
em | olB msec | msec 
177.38 1.89 +£0.05 30.3+0.5 87+ 2 
182.48 1.31 40.03 38.1+0.7 87+ 2 
187.27 0.78 +£0.02 50.2+1.2 sot 2 
189.22 0.566 + 0.016 

191.22 0.353 + 0.012 63 +4 85+ 6 
192.73 0.192+0.011 65 +8 78411 
193.76 0.105 + 0.010 

195 0 


Up to the time for these measurements the reactor had only run at very low power, 
so the background of photo neutrons was very low. The observed background mainly 
arose from the spontaneous fission of uranium and was less than one per cent even 
during the off-time of the pulse. 

The heavy-water level was read on the gauge glass in the control room simul- 
taneously with the temperature of the gauge glass and of the uranium rods. Corrections 
were made for the differences in temperature between the measurements. 

The critical heavy-water level was determined in the following way: A control rod 
was put into the reactor, the water was pumped up above the expected level and the 
reactor was run up to a power of a few watts. The water level was then lowered and 
the power kept constant by means of the control rod. This was continued until the 
control rod was completely out of the reactor. By small adjustments of the water 
the critical level could then be determined as accurately as a reading could be made 
on the gauge glass, that is + 0.2 mm. 

When the reactivity is to be calculated from a measured curve like the one shown 
in Fig. 2, some uncertainty may arise depending upon how the area is divided into 
the two parts A, and A,. This has been avoided by using eq. (A12) of Appendix IT, 
where directly measured values can be used in the numerical computation in both 
numerator and denominator. The figures used for the delayed neutron yields and 
periods are taken from Hughes et al. [16] and Bendt and Scott [17], and photo neutron 
data from Johns and Sargent [18]. The total fraction of delayed and photo neutrons 
was then 0.78 per cent. All groups of delayed and photo neutrons were assumed to 
have an energy of 0.5 MeV. By using eq. (A6) of Appendix I corrections were made 
for the difference in age between the prompt and delayed neutrons. An age difference 
of 9 cm? was assumed and the correction amounted to about 1 per cent for the 
largest reactivity and about 8 per cent for the smallest. In Table 1 the 9/8 values 
obtained from the measurements are listed together with estimated uncertainties. 


The influence of higher harmonic modes 


When evaluating the reactivities, corrections also had to be applied for the presence 
of higher harmonic modes in the spatial neutron distribution. This has been done 
by expanding the distribution of source neutrons in harmonic components and 
using eq. (A 26) in Appendix III. Since the target is located on one side of the reactor 
centre and the detector on the opposite side, the most serious harmonic mode is 
the one with a single radial nodal plane through the cylinder axis. This mode has a 
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Table 2. 
Distance from + 
reactor centre o/B = a 
cm 
44 0.89 49.1 
68 0.83 46.1 
87 0.85 47.7 
111 0.92 50.1 
136 0.86 47.3 
each + 0.02 each +1.4 


negative amplitude at the detector position, which means that a positive correction 
has to be applied to the reactivity. As an example, it may be mentioned that for the 
measurement shown in Fig. 2 a correction of two per cent was calculated. 

We also tried to study the influence of the harmonic modes by reactivity measure- 
ments in five different positions of the detector from close to the centre of the tank 
to the middle of the reflector. The results are given in Table 2. They have not been 
corrected for the higher harmonics, but the influence of the detector on the reactivity 
has been taken into account. (This was investigated with the pile running at a few 
watts, and the corrections were of the order of one per cent.) 

The scattering in the five values is greater than would be expected statistically. 
The large difference between the first value and the two subsequent ones can largely 
be ascribed to the most pronounced harmonic mode. However, the harmonic modes 
have a greater influence here than in the ordinary reactivity measurements, because 
the pulse length for these measurements was only one tenth of the repetition time, 
a circumstance which favours the higher harmonics. 


Uncertainties in the reactivity measurements 


To illustrate the uncertainties involved in the measurements we consider the case 
illustrated in Fig. 2, which, according to Table 1, gave the result 0/8 = 1.31 + 0.03. 

The statistical errors from the points in the curve give rise to an error of 0.2 per 
cent in the numerator and 0.8 per cent in the denominator of the right-hand side 
of eq. (A12). The periods and relative yields of the delayed neutrons enter in a similar 
way into both the numerator and the denominator and will therefore have a very 
small effect. We assume 0.2 per cent for the periods and 0.2 per cent for the yield, 
which means that the short periods could be incorrect by 10 per cent. The assumption 
of equal ages for all delayed and photo neutrons can at most give rise to 0.2 per cent 
error. This results in a 0.9 per cent error in (0/6 +1) or 1.6 per cent in e/f. A 20 
per cent uncertainty assumed in the correction for different ages between the prompt 
and delayed neutrons gives an error of only 0.2 per cent. The uncertainty due to 
the harmonic modes is estimated to be as large as the correction, namely 2 per cent. 
The error caused by non-saturation of the delayed neutron activities before measure- 
ments commenced could at a maximum amount to 0.5 per cent. The D—D neutron 
sources and variations in the source strength during the measurements do not in- 
troduce more than 0.1 per cent each. A total error of 2.6 per cent is obtained for o/Bp. 

If the absolute value of the reactivity is calculated we have in addition the un- 
certainty in 6, which mustbe considered to be as large as 6 per cent. 
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10 


Fig. 4. The reactivity asafunction Q5 
of 1/H? — 1/H2. The straight line 


was fitted to the measured points 
by the method of least squares. (6) 


(6) Of 02 03 04 05 06-1075 ¢m-? 
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The reactivity equivalent of the heavy water 


According to the one-group theory for a bare reactor we can write 


(6) 


ee — Wete a” M? ] 1 
CO kee 1+ MB\t 


where H is the effective height of the reactor and the index c denotes the value at 
criticality. 

In Fig. 4 the 9/8 values obtained from the measurements are plotted against 
1/H*? — 1/H:?. The effective critical height, H, was assumed to be 195 em according 
to the measurements by Larsson [3], which gave H,=195*}°. Theoretical work 
has given 200 em [3]. 

The points in Fig. 4 fit a straight line very well, and from the line we find the 
reactivity equivalent of the heavy water to be (7.2 +0.1) x 10-> per mm D,O at 
the critical level. This is in agreement with supercritical measurements by G. Jons- 
son [3], who obtained (7.0 + 0.1) x 10->. The errors do not include the 6 per cent 
uncertainty in the total yield of delayed and photo neutrons. 

From the slope of the line we find M? = 336 cm?, which at first seems too large. 
For the reactor lattice we can calculate M? = 235 from the measurements by Cohen 
[19]. According to calculations by Pershagen [3], the presence of reactor channels 
increases this by about 8.6 per cent. The presence of the reflector, not taken into 
account in eq. (6), raises the effective M? value by a further 25 per cent approximately 
[20]. We thus arrive at about 320 cm?. The remaining discrepancy can easily be 
explained by the uncertainty in the effective height of the reactor. While this does 
not have any influence on the reactivity equivalent of the water, an error of one 
per cent in the effective height will lead to as much as three per cent error in M?. 
By assuming the critical height to be 192 cm the M? value can thus be reduced to 
about 320 em?. The critical height is then still in good agreement with the experi- 
mental value. 
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Calibration of the safety rods 


Measurements similar to those described above were made in order to determine 
the reactivity equivalents of the safety rods. The heavy-water level was then 1.25 cm 
below the critical level, and measurements were carried out with safety rod number 
one, number two, or both rods in the reactor. The background was relatively higher 
for these measurements, as the reactor was farther from criticality, but it never 
amounted to more than 5 per cent. With no corrections for the harmonic modes we 
find after subtraction of the reactivity corresponding to the 1.25 cm of heavy water: 


Safety rod o/B @ (per cent) 
1 2.02 + 0.10 1.58 + 0.08 
2 2.06 + 0.10 1.614 0.08 
1+2 3.9 +0.3 3.1 +0.2 


As before, the uncertainty in f is not taken into account in the errors given. 

On another occasion, when the uranium lattice in the reactor was not quite 
complete, a few more measurements were made with the detector in two different 
positions. In the “‘in’”’ position, which was the same as that used previously, the 
detector was only about 30 cm from safety rod number one, while in all other cases 
the distance from the detector to the safety rods was at least 60 cm. The results are: 


o/B with the detector 
Safety rod 


in out 
1 2.22+ 0.06 2.00+ 0.07 
2 1.99+ 0.06 2.02 + 0.07 


These results do not indicate any strong dependence on the detector position. 
Because the neutron distribution is much disturbed by the safety rods it is difficult 
to calculate the harmonic corrections in this case. The harmonic mode having one 
single radial nodal plane through the cylinder axis will still be the most serious one, 
however. This means that a positive correction has to be applied to the reactivities. 
We assume 5 per cent for each safety rod and 10 per cent for both together. The 
final values are then: 


Safety rod Reactivity (per cent) 
1 1.7+0.2 
2 1.7+0.2 
1+2 3.44 0.4 


The results are in good agreement with the measurements by G. Jonsson [3]. 


Determination of the neutron lifetime 


In order to determine the mean lifetime of the thermal neutrons in the reactor a 
series of measurements was made with the pulse length only one tenth of the repeti- 
tion time, which was 0.25 or 0.5 sec. With this arrangement it was possible to follow 
the decay of the neutron density much more closely. These experiments were done 
at five of the heavy-water levels used previously. 
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As the repetition time was so long that the delayed neutron precursors decayed 
considerably during the off-time of the pulse, the measured curves had to be fitted 
to a sum of two exponentials. One corresponded to the prompt neutron decay, and 
one to that of the delayed neutrons. The fitting was done with a least-squares method, 
and the results are given in Table 1. 

The time constant /, of the prompt neutrons is related to the mean lifetime J of 
the thermal neutrons in the reactor according to eq. (5). By using this relation and 
the measured ratio o/f the values for //8 also given in Table 1 have been calculated. 
A small correction, less than 0.5 per cent, has been applied to eliminate the effect 
of the different sizes of the reactor during the different measurements. The values 
given are the mean lifetimes reduced to the critical size. 

The measurements described above, with the detector in different positions, were 
also used to study the influence of the higher harmonic modes on the time constant 
of the prompt neutrons. The results can be seen in Table 2. It is interesting to note 
that there is no tendency for the lifetime to increase as the detector is moved from the 
reactor centre to the reflector, a tendency of which Hellstrand [8] found a vague 
indication. In fact the spread between the values is not inconsistent with the 3 per 
cent statistical uncertainty in each measurement. It therefore seems as if the higher 
harmonic modes have a very small influence, which is also confirmed by the fact 
that the measurements of the mean lifetime made at different reactivities agree 
within the limits of error. 

As a weighted mean we get for 1/8 a value of 87 msec. Assuming 8 = 0.78 x 10-* 
we find a mean lifetime of 0.68 + 0.02 msec. The error here takes account of the un- 
certainty due to the harmonic modes but not of the 6 per cent uncertainty in f. 
The result agrees very well with the 0.71 + 0.03 msec obtained by Hellstrand [8], 
who used the slightly different 6 = 0.79 x 10-*. 

In this connection it may be pointed out how desirable it would be to have more 
accurate determinations of the total yield of the delayed and photo neutrons, f. 
The large uncertainty in this parameter enters into the absolute values of all the 
measured quantities in this paper. 


Conclusions 


The present investigation shows that the pulsed neutron source is a powerful 
means of studying a subcritical reactor. In addition to the measurements described 
above, the method can of course also be used for determination of the critical size 
of a subcritical assembly. As the uncertainties in the reactivity measurements de- 
scribed above are mainly of systematic type, it seems possible to increase the accuracy 
considerably. The most important thing would be to suppress the higher harmonic 
modes by a suitable arrangement of the source and the detector. With this improve- 
ment it would also be possible to make the measurements farther from criticality. 
It is therefore considered that the pulsed neutron source method applied to a sub- 
critical assembly may be a valuable complement to exponential experiments. 
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AppEnDIx I 


We consider a bare subcritical reactor which contains a neutron source. The source 
strength is assumed to be a periodic function of time, S(¢). It is further supposed 
that the reactor is so close to criticality that only the fundamental mode of the 
spatial neutron distribution is important. According to ref. [11], p. 295 ff. and with 
the same notation we can then write: 


Oe Lb AG C; —B%r 
renee se LYS AT CRE POT ERD 
ae fl—(1 Pykenlo+ > a Bry) * 5,04 BL) (Al) 
a0, k 
at = Aa Cait Gib Baa (Aa 


The ages of the prompt neutrons, the source neutrons, and the ith group of delayed 
neutrons are called t,, t;, and 7; respectively. As S(t) is a periodic function, ¢ and 
C;, must also be periodic functions when equilibrium has been reached for the delayed 
neutrons. When eqs. (Al) and (A2) are integrated over one period 7’, the terms on 
the left side therefore disappear and we obtain 


T T 7 
ats on = P — Bx f — Br, 
[1-—(l A) kaul | gat seas [22 | acide | sax), (A3) 
0 0 0 


T T 
k 
| accae~ “5.2. [ gat. (A4) 
0 3 0 
By summing all the groups of delayed neutrons we get from eq. (A4): 


T T 
Been | 40,40 *x, | pdt dpe ™, (A5) 
0 0 


T 
Now [odt can be eliminated between eqs. (A3) and (A5). The result is: 
0 


T 
— Br, 27, — Br: 
i f ese 1 — (1 — B) kote — here” P 2 Bie eis 


A yn if Kor = DB, eB ti ‘ (A6) 
ae io ef 1,0,dt _ i Pre 
0 


i 


In the case when all groups of the delayed neutrons have the same age as 
the prompt neutrons, this expression reduces to 
A, ay 1 4 kere Q 

A, *phageed (A7) = (3) 


242 


ARKIV FOR FYSIK. Bd 11 nr 13 
ApPpENDIx II 


If eq. (A2) is multiplied by exp (4,t) and integrated over the period 7, 
we obtain 


T 
CT) eT — 00) =F 8,2 | ged, (A8) 
ei 0 
As C; is a periodic function, C;(7')=C;(0), and thus 
ake 
Oy Ss i At 
C; (0) AE zis, | pet’ dt. (A9) 


Multiplying by exp (—B?r;) 4; and summing over the delayed neutron groups, 
we get 


L627 4t1,0, (0) = ies fei a 7 vee (A 10) 


i p 2 


Combination with eq. (A5) gives 


ne *% 2,0; (0) | pdtd pe” 


v 


zern | acdt~ : (All) 
i eo Bt 
UB A oH ae | ve dt 


By substitution in eq. (A3) we can obtain a corresponding expression for 
qr 


exp (—B?x,) f Sdt. The ratio is 
0 


T 


ern Sdt 


A, 0 


ze 
a | A,C, dt 
0 (A 12) 


T 
ert at 
UB AF pe dt 
= a —1. 


fn — Bt NBS ee, 
Pe(acliekdeei1 ae 
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Appenpix III 


If the treatment in chapter 7 of ref. [11] is extended to include one group 
of delayed neutrons, we get the following equations 


; 29) § Dy*B (7,1) Lah (1) + PI Tet) (A 13) 
viq(r, 1.) = “EE (A14) 

éC(r, 0 k 
Ane RU eta Aa (A 15) 
4 (7, 0,1) = Ze(1—B) p(r, N+ 2C(% +8 (0,0. (A 16) 


We have here assumed that the prompt, delayed, and source neutrons have 
the same energy. By separating the variables in eq. (A14) and introducing 
the eigenfunctions obtained from the boundary conditions we get 


(r= 3 Ky Bale Be Tal (A17) 

S(r,)= > Sa) Ral) (A118) 

C(r, t) -> C,, (t) Ry (7) (A 19) 

b= ps tm S KaPal)—ACo()— Sa} Bal). (20) 


These expansions are put into eqs. (A13) and (A15). Integrating over one 
period, 7’, and simplifying yields the following relations: 


T T ; 
J S,dt=(1-k,) K,f T,dt, (A21) 
0 0 

T Wh 

JAC, dt=Bk,K,{ Tydt. (A 22) 
0 0 


By sptea of eq. (A20) over one period, and using the above expressions, 
we ge 
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T 
xe t)dt= Ss Det hy m| 4 iia) (A 23) 
0 


If now eqs. (A21) and (A22) are added and K,, r T,,dt is eliminated, we find 
0 


the following expression for the area under a measured curve: 


T T 
| eeoae- = 2 tial {+ dt+ | s ai R( (r).  (A24) 
0 


0 


The ratio of the two parts of the area is then 


T 
foe] kn 
A Page aay) B) oJ S,,dt 
1 
zi 4 (A 25) 
o In 
Somat | Acnae 
0 
This is finally transformed to 
T 
* hk, 1k (I ean: 
i n\T) 0 
Lo 
ay 1-k (Lp) R, (r) T 
[S,dt 
2 eg 
As - 4 (A 26) 


4 

en se IZ, 1—-& ey J 
ght, 1— kyl =, (1—F} Pay 

0 
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